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Abstract 


An  iterative  solution  to  the  problem  of  scattering  from  a  one-dimensional  rough  surface 
is  obtained  for  the  Dirichlet  boundary  condition.  The  advantages  of  this  method  are  that 
bounds  for  the  convergence  of  the  solution  may  be  established  and  that  the  solution  may  be 
readily  iterated  to  sufficiently  high  order  in  the  interaction  to  examine  the  rate  at  which  it 
converges.  Absolute  convergence  of  the  iterative  solution  is  also  a  sufficient  condition  for  the 
convergence  of  the  operator  expansion  method  for  surfaces  on  which  the  slope  is  everywhere 
less  than  unity.  A  numerical  example  of  scattering  from  an  echelette  grating  is  considered, 
and  bounds  for  convergence  established.  It  is  found  that  for  scattering  from  such  surfaces, 
the  rate  at  which  the  iterative  solution  converges  decreases  as  the  surface  slope  is  increased. 
Corresponding  results  are  found  for  the  operator  expansion  method. 
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1.  Introduction 

In  the  past  two  decades,  many  new  approaches  to  the  rough-surface  scattering  problem 
have  been  proposed  that  have  potential  for  treating  scatter  from  very  rough  surfaces.  These 
include  the  small-slope  approximation  (Voronovich  1985)  and  its  extensions  (Tatarskii  1993, 
McDaniel  1995,  and  Charnotskii  and  Tatarskii  1996).  Also  noteworthy  in  this  respect  are 
the  phase  perturbation  technique  of  Shen  amd  Maradudin  (1980)  and  the  development  of 
second  order  corrections  to  the  Kirchhoff  approximation  by  Chen  and  Ishimaru  (1990).  The 
operator  expansion  method  of  Milder  (1991)  may  also  be  applicable  to  very  rough  surfaces. 
This  method  is  basically  an  extension  of  earlier  work  by  Lysanov  (1956)  and  Meecham 
(1956)  to  obtain  higher  order  terms  of  a  series  solution  to  the  surface  scattering  problem 
(Milder  1996).  In  a  further  extension,  Voronovich  (1996)  has  exploited  the  fact  that  the 
operator  expansion  method  is  not  gauge  invariant  to  formulate  the  non-local  small-slope 
approximation. 

It  is  not  surprising  that  operator  expansion  methods  are  more  accurate  than  small-slope 
approximations  since  the  small-slope  approximations  are  obtained  through  the  expansion  of 
certain  terms  in  the  operator  expansions  (McDaniel  1995).  For  scatter  from  random-rough 
surfaces,  the  small-slope  approximations  permit  the  calculation  of  ensemble  averages,  while 
the  operator  expansion  methods  require  Monte  Carlo  simulation.  However,  such  simula¬ 
tions  may  be  readily  performed  through  the  repeated  application  of  the  Fourier  transform 
(Kaczkowski  and  Thorsos  1994).  An  important  remaining  issue  is  the  regime  of  apphcability 
of  operator  expansions. 

In  this  article,  an  iterative  solution  to  an  integral  equation  of  the  first  kind  for  the  normal 
derivative  of  the  surface  field  is  obtained  for  scattering  from  a  one-dimensional  randomly 
rough  surface  characterized  by  the  Dirichlet  boundary  condition.  This  results  in  a  series  for 
the  surface  potential  for  which  the  conditions  for  convergence  may  readily  be  established. 
Successive  terms  in  this  series  may  also  be  easily  generated  to  examine  the  rate  at  which 
the  series  converges.  A  disadvantage  of  this  solution  is  that  obtaining  numerical  results 
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requires  even  more  computational  effort  than  in  the  case  of  the  operator  expansion  method. 
In  addition  to  the  need  for  Monte  Carlo  simulation  for  treating  randomly  rough  surfaces, 
double  integrations  are  required  for  scatter  from  one-dimensional  surfaces,  and  four-fold 
integrations  for  the  two-dimensional  case.  Because  operator  expansions  are  obtained  by 
series  expansion  of  terms  in  the  iterative  solution,  a  sufficient  condition  for  the  convergence 
of  the  operator  expansion  method  may  be  obtained  for  surfaces  having  slopes  everywhere 
less  than  unity  where  this  series  expansion  is  valid. 

In  section  2  of  this  paper,  the  basic  theoretical  development  of  the  iterative  solution  is 
presented  and  its  relationship  to  the  corresponding  operator  expansion  is  discussed.  The 
surface  potential  for  a  plane  wave  incident  on  a  flat  surface  of  arbitrary  slope  is  evaluated  to 
demonstrate  the  differences  between  the  iterative  solution  and  operator  expansion.  Section 
3  addresses  the  numerical  solution  of  the  relevant  equations  and  establishes  the  conditions 
for  convergence  of  the  iterative  and  operator  series. 

Scattering  from  a  echelette  grating  is  considered  in  section  4,  in  which  numerical  results 
for  the  iterative  solution  and  operator  expansion  are  compared  with  exact  solutions.  First, 
the  region  of  convergence  of  the  iterative  solution  is  established  and  is  found  to  significantly 
exceed  that  of  Milder’s  (1991)  Fresnel  phase  criterion.  The  rate  of  convergence  of  the  iterative 
solution  is  examined  by  computing  successive  terms  of  this  series  and  is  found  to  depend  on 
the  incident  grazing  angle  as  well  as  the  surface  slope.  Numerical  results  are  also  presented 
for  the  operator  expansion  method.  In  section  5,  the  findings  of  this  study  are  discussed. 

2.  Theoretical  development 

The  problem  addressed  is  illustrated  in  figure  1  in  which  a  plane  wave  of  wavelength 
A  =  2t:  jk 

Pine  =  exp[?(«oa:  -  i/qz)]  (1) 

is  incident  on  a  one- dimensional  rough  surface  on  which  the  Dirichlet  boundary  condition  is 
obeyed.  Here  Kq  =  kcosOo,  I'o  =  fcsin^o?  is  the  incident  grazing  angle,  and  k  is  the  wave 
number. 


2 


/dx 


(2) 


(3) 


The  scattered  field  p,  is  represented  as 

Pa  =  J  dKfS{Kf,Ko)exp[i{KfX  +  l/fZ)] 

where  Kf  =  kcosOf,  vj  =  ksinOf,  and  the  scattering  amplitude  S{Kf,Ko)  is  defined  as 

dx 

where  the  surface  potential  ^{x)  is  defined  below.  The  definition  of  S{Kf,Ko)  in  (3)  differs 
from  that  of  Voronovich  (1985),  and  the  requirement  that  it  be  reciprocal  takes  the  form 

UfS{Kf,Ko)  =  i/o«S'(-Ko,-«/). 

Evaluation  of  the  scattering  amplitude  requires  the  solution  of  an  integral  equation  (Uret- 
sky  1965)  for  the  surface  potential 


J  il){x')V{x,x')dx'  =  —2tp[x) 


where 


tp{x)  =  - 


dp 


dp 


Z=C(®) 


V{x,x')  =  kH^o’{kp)l2 


(4) 


(5) 

(6) 


with  p  =  {(a;  — a;')^  +  [C(a;)  — C(a:')]^}^/^,  a  Hankel  function  of  the  first  kind,  ('  the  surface 

slope,  and  ^(x)  =  pinc{x).  In  (4),  ^(x)  is  evaluated  on  the  rough  surface  ^(x).  Applying  the 
operator  f  dx  exp(—iKx)  to  both  sides  of  (4)  yields 
dxdx' 


/dxdx^  t  ^ 

— 2^V(x,x')exp(— i«x)  j  dK'tp(K')exp{iK'x')  =  — 2^(/c) 

where  V’(«)  =  f  dxtl>{x)exp{—iKx),  and  a  similar  expression  holds  for  ^(«). 
To  proceed,  we  write 


(7) 


(8) 


V(x,x')  =  -  V<')(x,x') 

where  V(“)(x,x')  =  kH^^\k\x  -  x'|)/2,  and  VW(x,x')  =  V(°)(x,x')  -  V(x,x').  The  surface 
potential  V’(«)  is  next  expanded  as  an  ordered  series  in  the  interaction 

V’(«)  =  ^  +  •  •  •  •  (9) 
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Substituting  (8)  and  (9)  in  (7),  and  solving  successively  for  the  yields 


j  ^^^^(/c')A(«,  K')d«'  =  — 2^(/c) 

(10) 

where 

A(k,  k)  =  f  0®^(x,  X  )  exp[i(/cV  —  Kx)]. 

J  ZTT 

(11) 

The  integral  on 

the  left  hand  side  of  (10)  may  be  evaluated  to  obtain 

^(^1(«)  =  —2u^{K)fk 

(12) 

where 


Higher  order  terms  for  which  n  >  1 ,  take  the  form 

k^^'^\K)lv  =  j  K')dK'  =  I  t/>("-^)(/c')5(/C,  K')dK'  (13) 

where  B{k,k')  is  defined  similarly  to  A{k,k')  in  terms  of  Uretsky  (1965)  has 

shown  that  multiple  integrals  of  the  type  appearing  on  the  left  hand  side  of  (7)  converge 
provided  V’(«0  grows  no  faster  than  some  finite  power  of  so  that  the  interchange  of  the 
integral  over  k'  with  the  spatial  integrals  to  obtain  (10)  and  (13)  is  permissable. 

To  obtain  the  operator  expansion,  -  H^^\k\x  -  x'\)  is  expanded  using  (Watson 

1966) 

i/W|z(l  +  r)V^|  -  =  E  ^  (^)”  (11) 

The  convergence  of  this  series  is  governed  by  the  singular  nature  of  the  Hankel  function. 
Inserting  the  limiting  form  of  the  Hankel  function  into  (14),  it  follows  that  the  series  is 
absolutely  and  uniformly  convergent  provided  r  =  [(^(x)  -  C(a:')]^/(®  -  x'Y  <  1,  or  that,  the 
surface  slope  is  everywhere  less  than  unity.  Substituting  (14)  into  (9)  yields  an  expansion 
for  the  surface  potential.  If  (9)  and  (14)  are  absolutely  convergent,  the  summation  of  this 
expansion  may  be  rearranged  in  any  arbitrary  order  since  the  orders  of  integration  and 
summation  in  (9)  may  be  freely  interchanged.  In  particular,  the  sum  may  be  rearranged 
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to  form  a  simple  series  in  powers  of  the  surface  excursion  to  obtain  an  operator  expansion. 
Thus,  under  these  conditions,  the  operator  expansion  will  converge  absolutely. 

The  operator  expansion  may  also  be  obtained  by  substituting  (14)  directly  into  (7). 
Expanding  the  surface  potential  in  an  ordered  series  in  the  surface  elevation 

V’(«)  =  '0l(«)+V’2(«)  +  V'3(«)  +  ---  (15) 

and  also  inserting  it  into  (7),  yields  an  expression  that  may  be  solved  successively  to  de¬ 
termine  the  i/’t-  In  this  case,  the  spatial  integrals  can  be  performed  to  obtain  (McDaniel 
1995) 

■0i(/c)  =  —2v^[K,)lk  (16) 

and 

2V>2(«)  =  y  dK'xpi{K')6{K'  +  q  +  q'  -  k)  j  dqdq'C{q)C{q')F{q,  q')  (17) 

where 

+  <l  +  <l')  +  *^('«0  ~  +  9)  -  +  q') 

and  ({q)  is  the  Fourier  transform  of  the  surface  displacement  C(x) 

C(9)  =  /  ^C(®)exp(-iga:). 

On  comparing  (12)  and  (16),  it  is  evident  that  the  leading  terms  of  the  iterative  solution 
and  the  operator  expansion  are  identical. 

Because  F{0,q')  =  F{q,0)  =  0  in  (17),  it  may  be  argued  (Voronovich  1985)  that  V’2(k)  is 
of  second  and  higher  order  in  the  surface  slope.  Thus,  following  Voronovich  (1985,  1996),  an 
integration  by  parts  may  be  performed  to  extract  the  contributions  of  second  order  in  the 
slope  from  (17) 

2^2('«)  =  ik  J  dK' —  k'){v  +  1/0  —  v{k')  —  v{k  —  /c')]/i/(/c').  (18) 

This  result  will  be  used  in  examining  the  convergence  of  the  operator  expansion  method. 
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The  differences  between  the  iterative  and  operator  series  are  apparent  even  in  very  simple 
cases  such  as  scattering  from  a  flat  surface  of  slope  s  for  which  C(x)  =  sx.  The  exact  solution 
for  the  surface  potential  is 

i/;e(x)  =  — 2exp(i«x)(P  —  jk  (19) 

where  P  =  +5^)  and  k  =  Kq  —  t'O'S-  The  iterative  solution  yields 

=  — 2exp(iKx)(P  —  P)^^^[2  —  {k^  —  k^Y^^f(k^  —  k^Y^^]/k.  (20) 

Finally,  the  operator  expansion  yields 

^i(x)  +  =  —  exp(iKx)(P  —  k^)^/^[2  +  s^k^l{k^  —  k^)]/k.  (21) 

The  surface  potential  for  the  iterative  solution  (20)  may  be  obtained  from  (19)  by  writing 
(P  —  P)^/^  =  (P  —  P)^/^  +  (P  —  P)^/^[l  —  (P  —  P)^/^/(P  —  P)^/^]  and  expanding  the 
coefficient  of  the  second  term  on  the  right  hand  side  of  this  expression  for  small  s.  Further 
expansion  of  (20)  for  small  s  yields  (21).  Because  the  next  higher  order  terms  in  these 
expansions  are  of  fourth  order  in  the  slope,  both  (20)  and  (21)  are  valid  to  third  order  in  the 
slope. 

The  behavior  of  (20)  and  (21)  as  a  function  of  the  incident  grazing  angle  $q  differs 
significantly:  (21)  is  unbounded  when  Oq  or  Of  vanishes  irrespective  of  the  slope,  while  (20) 
vanishes  at  these  points.  On  the  other  hand,  (20)  is  unbounded  for  glancing  incidence 
on  the  surface  where  tan^o  =  and  ipE  =  0.  Because  higher  order  terms  in  the  iterative 
solution  contain  successively  higher  order  powers  of  (P  —  /c^)~^/^,  it  is  clear  that  the  iterative 
solution  cannot  be  truncated  at  arbitrary  order,  for  this  example  where  the  surface  excursion 
is  unbounded.  It  is  thus  anticipated  that,  not  only  the  slope,  but  the  rough  surface  excursion 
as  well,  will  play  a  role  in  the  series  convergence. 

3.  Numerical  solution 

This  section  addresses  the  numerical  evaluation  of  the  iterative  solution  and  the  con¬ 
ditions  for  its  convergence.  The  random-rough  surface  ^(a:)  is  treated  as  periodic  with  a 
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fundamental  period  A  so  that  scatter  occurs  at  grating  angles  0^  where 


cos  On  =  cos  ^0  +  27rn/(/;A). 


(22) 


It  is  convenient  also  to  define  =  k  cos  On  and  t/n  = 

The  integral  equation  for  the  surface  potential  (7)  is  replaced  by  the  matrix  expression 
(Holford  1981) 

=  — 2'0 


where 


^n=j^  ^PinceXp{-iKnX) 


A 
dx 


/*A  /*oo 

Mmn=  -T-  dx'V{x,x')cX.p[i{KnX' -  KmX)]. 
JO  A  J^-oo 


Using  (8)  to  define  A  and  B  yields 


^mn  kSfn—nl^v 


B  =  A  -  M. 


With  (24)  and  (25),  (23)  may  be  rewritten  as 


(23) 


(24) 

(25) 


•0  —  A  ^BV»  =  — 2A 


(26) 


The  exact  solution  to  (26)  is 


= -2(1  -  C)-^A-'V’ 


(27) 


and  the  iterative  solution  is 


ij}  =  -2(1  +  C  +  +  •  •  OA-^^ 


(28) 


where  C  =  A  ^B. 
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The  conditions  under  which  rp ^  converges  as  the  dimensions  of  the  matrix  C  are  increased 
have  been  investigated  by  Urusovskii  (1965)  and  McDaniel  and  Krauss  (1991).  The  concern 

here  is  the  condition  under  which  the  series  1  +  C+C^+C^H - converges  to  (1  —  C)“^.  The 

requirement  for  convergence  is  that  the  spectral  radius  of  C  be  less  than  unity  (Varga  1962), 
where  the  spectral  radius  p{C)  =  |Ai|  and  Aj  is  the  eigenvalue  of  the  n  X  n  complex  matrix 
C  of  greatest  absolute  value.  In  this  case,  the  series  ipn  =  converges  absolutely. 

Numerical  studies  are  necessary  to  further  pursue  conditions  for  the  convergence  of  the 
iterative  and  operator  series. 

4.  A  numerical  example 

Numerical  results  are  obtained  for  scattering  from  a  periodic  grating  described  by 


f  7ra(l  —  4x/A)/2;  0  <  x  <  A/2  facet  1 

I  7ra(— 3  +  4x/A)/2;  A/2  <  x  <  A  facet  2. 


(29) 


The  slope  of  the  grating’s  facets  tane  =  27ra/A,  and  its  rms  roughness  cr  =  Traf\/l2. 

First,  let  us  consider  the  range  of  the  parameters  kcr  and  tane  for  which  the  iterative 
solution  converges.  For  this  surface,  the  elements  of  the  matrix  M  derived  in  the  appendix 
are  given  by 


/_2x  ai/n+j[l  -  (-1)^'  exp(ii^n+jair)] 

■Mmn  —  {^ka/ir  )  2^  j—  —  I  — -  — r  - 


(30) 


Figure  2  shows  values  of  the  two  parameters  for  which  p{C)  is  unity  for  various  angles  of 
incidence.  For  tan  e  <  1.5,  aU  of  the  points  lie  to  the  right  of  the  line  kcr  =  irf  tan  e  suggesting 
that  this  line  is  an  upper  bound  on  the  parameters  that  will  always  yield  convergent  iterative 
solutions  for  low  to  moderate  surface  slopes. 

Figure  3  shows  values  of  the  spectral  radius  computed  along  the  line  fccr  =  tt/  tan  e  for  two 
values  of  the  incident  grazing  angle  9o  =  7r/2,  7r/4.  Here,  again  it  appears  that  for  surfaces 
of  low  to  moderate  slope  a  sufficient  condition  for  convergence  of  the  iterative  solution  is 
that  A:crtane  <  tt. 


8 


For  slopes  of  less  than  unity,  the  operator  expansion  method  will  also  converge  in  this  re¬ 
gion  by  virtue  of  the  arguments  presented  in  sections  2  and  3.  This  bound  may  be  applicable 
to  scattering  from  other  surfaces  as  well.  Kaczkowski  (1994)  applied  the  operator  expansion 
method  to  random-rough  surfaces  characterized  by  a  Gaussian  wavenumber  spectrum.  For 
a  surface  having  an  rms  slope  a'  =  0.85,  and  an  incident  grazing  angle  of  60°  he  found  for 
k(Tcr'  =  3.56  that  the  method  converged  for  40°  >  ^/  >  140°  and  appeared  to  diverge  for 
other  scattering  angles.  Kaczkowski  also  reported  that  for  higher  values  of  kcrcr\  the  operator 
expansion  method  failed  to  converge. 

If  the  iterative  solution  or  the  operator  expansion  are  to  be  useful  in  practical  applications, 
they  must  not  only  converge,  but  converge  rapidly.  To  examine  the  convergence  of  these 
series,  an  example  presented  by  Voronovich  (1996)  is  considered  for  which  tane  =  0.318  and 
^0  =  0.609.  In  this  case,  ka  has  been  reduced  to  9.4  to  insure  convergence.  The  reflection 
coefficients  shown  in  figure  4  obey 

-  ^  rA 

Rn  =  2xS{Kn,Ko)  =  {k/2i/n)Y^ipm  exp[?x(/c„i  -  «n)  -  ii/„C(a:)].  (31) 

m  Jo  A 

In  figure  4,  predictions  of  the  iterative  solution  are  compared  to  the  exact  solution  -  with  the 
numbers  in  parentheses  denoting  the  order  in  the  slope  to  which  the  predictions  are  valid. 
It  is  evident  from  this  comparison,  that  terms  of  5th  order  in  the  slope  must  be  retained  to 
obtain  good  agreement  with  the  exact  solution. 

Corresponding  results  are  shown  for  the  operator  expansion  in  figure  5,  in  which  the 
numbers  in  parentheses  again  indicate  the  order  in  the  slope  of  the  predictions.  For  the 
periodic  surface  considered  here  (17)  and  (18)  take  the  form 

^  CnCm— n— ?' (^m  “h  ^j+n 

j,n 

j 

where 

dx 

Cn=  —Cix)  exp{—2irinx  /  A)  =  2a/{irn^); 

J  0 


n  odd. 


The  operator  expansion  appears  to  converge  more  rapidly  than  the  iterative  solution  near 
the  forward  scattering  peak  at  n  =  0  which  corresponds  to  double  scattering,  first  by  facet 
one  and  thence  by  facet  two.  The  peak  at  n  =  -15  is  due  to  single  scatter  from  facet  two. 
Neither  series  converges  rapidly  when  the  scattered  field  is  near  grazing  where  n  >  0. 

The  predictions  of  the  iterative  solution  for  scattering  from  a  surface  on  which  tan  e  is 
greater  than  unity  are  shown  in  figure  6.  In  this  example,  ka  —  2  and  tane  =  1.2.  The 
magnitude  of  the  reflection  coefficient  R_i  converges  rapidly  when  the  incident  field  is  near 
grazing,  Bq  <  30°,  but  at  higher  angles,  terms  of  ninth  order  in  the  slope  must  be  retained 
for  agreement  with  the  exact  solution.  For  completeness,  the  corresponding  predictions  of 
the  operator  expansion  method  are  shown  in  figure  7.  In  this  case  there  is  no  gaurantee 
that  results  are  convergent,  however,  retaining  terms  of  second  and  third  order  in  the  slope 
appears  to  provide  better  agreement  with  the  exact  solution  -  especially  at  higher  grazing 
angles. 

5.  Discussion 

This  paper  has  addressed  the  two-dimensional  problem  of  scattering  from  a  rough  surface 
on  which  the  Dirichlet  boundary  condition  is  imposed.  The  method  of  successive  approxima¬ 
tions  has  been  applied  to  an  integral  equation  of  the  first  kind  to  obtain  an  iterative  solution 
for  the  surface  potential.  The  conditions  under  which  this  iterative  solution  converges  have 
been  established,  and  it  has  been  shown  that  this  regime  of  convergence  is  shared  also  by 
the  operator  expansion  method  with  the  only  additional  provision  being  that  the  slope  of 
the  rough  surface  is  always  less  than  unity.  While  the  conditions  for  convergence  obtained 
in  a  numerical  example  are  generous  compared  to  Milder ’s  (1991)  Fresnel  phase  criterion 
k(T(j'  <C  1,  where  a'  is  the  rms  surface  slope,  they  are  still  quite  restrictive.  Milder ’s  criterion 
may,  however,  play  a  role  in  the  rate  at  which  operator  expansions  converge. 

McDaniel  (1995)  has  shown  that  small-slope  expansions  may  be  obtained  by  expanding 
terms  of  the  type  exp{ii/o,/[C(ic)  —  C(®0]}  operator  expansion’s  scattering  amplitude. 

The  arguments  of  section  2  are  not  applicable  in  this  case,  and  further  studies  are  necessary  to 
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determine  the  range  of  parameters  for  which  this  series  converges.  The  last  approximation 
that  we  consider  in  this  hierarchy  of  solutions  is  perturbation  theory,  which  results  from 
expanding  the  term  exp[i{uf  +  i^o)C(aj)]  of  the  scattering  amplitude  obtained  from  of  the 
small-slope  approximation.  The  arguments  of  section  2  are  again  inapplicable,  however,  it  is 
clear  in  this  instance,  that  the  rate  of  convergence  of  the  solution  will  be  strongly  dependent 
on  the  magnitude  of  the  exponential’s  argument. 

Currently,  the  region  of  convergence  of  these  last  two  solutions  and  the  rate  at  which  they 
converge  can  only  be  determined  through  numerical  studies,  which  are  severely  limited  by 
the  complexity  of  higher  order  terms  in  these  approximations.  For  example,  the  perturbation 
based  small-slope  approximation  requires  the  evaluation  of  27  terms  to  obtain  a  result  valid 
to  third  order  in  the  slope  (Thorsos  and  Broschat  1995)  and  McDaniel’s  (1995)  extended 
small- slope  approximation  contains  more  than  100  terms  at  this  order  in  the  slope.  Because 
higher  order  terms  in  the  iterative  solution  may  be  readily  generated,  numerical  studies  based 
on  this  series  may  provide  insight  into  the  convergence  rates  of  the  other  approximations. 
As  an  example,  the  slow  convergence  in  figure  4  for  n  >  0,  may  lead  to  slow  convergence  of 
the  operator  and  smaU-slope  series  when  the  scattered  field  is  near  grazing. 
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Appendix 

The  elements  of  the  matrix  M  are  evaluated  as  in  Uretsky’s  (1965)  treatment  of  scattering 
by  a  sinusoidal  surface.  Here,  the  Hankel  function  appearing  in  the  matrix  will  be  represented 
by 

/-v  too 

H^\kp)=  j  —  exp[ffc:^(a:-a;')-t-ij/^|C(a;) -C(«')l]  (Al) 

J—oo  irUx 

where  Noting  that  exp[fi/i.|^(a;)  —  C(®0I]  ^  periodic  function  of  x'  for  fixed  x 

it  may  be  expressed  as 

exp[ii/^|C(x)  -  C(a;')l]  =  £  exp{2Trij  x' I  A)  (A2) 

3 
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where 


/•A 

Cj{x,  Ux)  =  —  exp[zi/^|C(x)  -  C(x")|  -  27rijx"/A]. 

The  matrix  elements  Mmn  then  take  the  form 

rA  (lx 

Mmn  =  {k/2x)  —exp{-iKmx)  dx' exp{iKnx') 

Jo  A  J’-oo 


/oo 

— -  exp[j/:j;(x  -  x')]  ^  Cj{x^  v^)  exp(27rijxYA). 

-oo 


dkx 

-oo  l>x  j 

Performing  the  integral  over  x'  in  this  expression  yields 

dx 


M„ 


tlx  d^x 

Jq  X  “  ”^)/A]C'j(x,  Un+j)ll'n+j 


(A3) 


(A4) 


(A5) 


where  the  interchange  of  the  summation  and  integration  to  obtain  (A5)  has  been  justified 
by  Uretsky  (1965).  Using  (A3)  in  (A5),  the  spatial  integrals  may  readily  be  performed  to 
obtain 

Mmn  =  M21+M!^2 


where 


=  {2ka/ir^)  77 


ai/n+j[l  -  (-l)^'exp(ii/n+ja7r)] 


f  \j^  -  «^t'^+j][(n  -m  +  jy  - 
= -{ikaM  E  ■ 

j  J  ^n+j 
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Figure  1 .  Geometry  for  scattering  from  a  random-rough  one-dimensional  surface. 
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Figure  2.  Values  of  kcr  and  tane  for  which  the  spectral  radius  /3(C)  is  unity. 
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Figure  4.  Magnitude  of  the  reflection  coefficient  as  a  function  of  reflected 
order  predicted  by  the  iterative  solution  for  kcr  =  9.4,  tane  =  .318,  and  an 
incident  grazing  angle  =  .609.  The  numbers  in  parentheses  indicate  the 
order  of  the  slope  to  which  the  predictions  are  valid. 
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Figure  5.  Magnitude  of  the  reflection  coefficient  as  a  function  of  reflected 
order  predicted  by  the  operator  expansion  for  the  case  of  figure  4. 
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Figure  6.  Magnitude  of  the  reflection  coefficient  i?_i  as  a  function  of  in¬ 
cident  grazing  angle  predicted  by  the  iterative  solution  for  ka  =  2  and 
tane  =  1.2.  Again,  the  numbers  in  parentheses  indicate  the  order  of  the 
slope. 
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Figure  7.  Magnitude  of  the  reflection  coefiicient  i?_i  as  a  function  of  inci¬ 
dent  grazing  angle  predicted  by  the  operator  expansion  for  the  case  of  figure  6. 
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